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Abstract. We present a detailed perturbative study of the dynamics of several types of atom- 
atom correlations in the famous Fermi problem. This is an archetypal model to study micro- 
causality in the quantum domain where two atoms, the first initially excited and the second 
prepared in its ground state, interact with the vacuum electromagnetic field. The excitation can 
be transferred to the second atom via a flying photon and various kinds of quantum correlations 
between the two are generated during this process. Among these, prominent examples are 
given by entanglement, quantum discord and nonlocal correlations. It is the aim of this paper 
to analyze the role of the light cone in the emergence of such correlations. 



1. Introduction 

Since its original conception, the two-atom Fermi problem [1J has been the subject of an 
intense academic debate. Stated in very simple terms, this gedanken experiment goes as 
follows: suppose there are two two-level atoms (qubits), say A and B, spatially separated by 
a distance r and interacting independently with a the (multi-mode) quantized electromagnetic 
field, initially in the vacuum state. The atom A is in an excited state |e) whereas the atom B 
is in its ground state \g). At some time t the atom A emits a photon. The original question 
by Fermi was then the following: as long as the two atoms are causally disconnected, is 
the excitation probability of B independent on the presence of A7 This question points at 
the very foundations of quantum mechanics: do quantum mechanical probabilities respect 
micro-causality? Over the past years a considerable amount of literature has dealt with this 
problem and several perturbative solutions have been proposed [|2l [3l HI [51 0, sometimes 
even presenting opposite conclusions. Recently, in Ref.[|7]|, a non-perturbative proof of 
strict causality in the Fermi problem has been finally given along with an explanation of 
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why the existence of correlations outside the light-cone connecting the two atoms is not in 
contrast with micro-causality. In fact, in a previous paper flS), some of the same authors had 
already studied the dynamics of concurrence [9J and found that it starts increasing just before 
(t — to) = r/c by a very tiny amount. 

In this paper, motivated by such results, we take a step further and investigate the 
dynamics of other types of atom-atom correlations. In particular, besides extending the 
analysis of entanglement dynamics, we study the time evolution of geometric quantum discord 
and (classical) connected correlation functions. Quantum discord was first introduced by 
Ollivier and Zurek [TTOl fTTTl as a novel measure of the quantumness of correlations. The idea 
behind it is conceptually very simple. Quantum discord is indeed defined as the discrepancy, 
in the quantum regime, between two classically equivalent definitions of mutual information. 
It is believed to capture a more general type of quantum correlation than entanglement, in 
the sense that quantum states with zero entanglement but non-zero quantum discord do exist 
(see e.g. lO^[T^[T4l[T5l[T6l[T7l '). Unfortunately, the computation of quantum discord implies 
solving a rather complicated minimization problem and, although considerable improvements 
have been made over the last years [fT8l . analytical results are available for very few cases 
only |fT9l |20l |2"U| . In order to overcome such computational issues, alternative indicators of 
general quantum correlations have been recently introduced, mostly based on distance-based 
approaches [|22ll23l . On the other hand, connected correlation functions provide a statistical 
quantifier of the (classical) correlations extractable during a joint measurement of the two 
atoms Il24l|25ll26t 

The results reported in this article have a double merit: on one hand they give a more 
complete picture of the dynamics of correlations in the two-atom Fermi problem. In fact, 
to the best of our knowledge, this is the first study of quantum discord and more general 
correlation dynamics in such a physical model, where an exact solution of the dynamics is 
still missing. As it turns out, all the types of correlations we consider have a nice physical 
interpretation in terms of a few relevant physical processes of the dynamics. On the other 
hand our results suggest a way to detect atom-atom correlations outside the light cone, that is, 
when the two atoms are causally disconnected. It is important to remark that our calculations 
are performed in the framework of time-dependent perturbation theory and they are exact and 
consistent up to the second order in the coupling constant. However, this is certainly not a 
problem since, as stated above, strict causality in the model has been analytically proven with 
no usage of perturbation theory [7J. Moreover, the time interval under scrutiny falls within 
the limits of validity of our approach provided the two atoms are not that distant. 

2. The Model 

We will consider a one-dimensional physical setup. A pair of two-level superconducting 
qubits (artificial atoms) A and B, separated by a fixed distance r, interact with an 
electromagnetic field propagating along the open transmission line connecting them. We 
name the atomic levels as {\g), |e)} and we assume the following multimode structure for 
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the field 

V(x) = J dk^[N^ k [e ikx a k + e" ifcl 4] (1) 

where iV is a normalization factor which may accommodate different circuit QED 
architectures, the dispersion relation is linear u k = v\k\, and a k , a\ are the usual annihilation 
and creation operator satisfying boson commutation relations [a k , a\.,] = 5 kjk '. We define 
Vtj = coj e — ujjg (J = A, B) the energy separations between the qubit levels and we assume 
the qubits to be much smaller than the relevant wavelengths Aj = v/(£lj/(2n)), v being the 
propagation velocity of the field quanta which in this scheme depends on the microscopic 
details of the model. Specifically, v — 1/ y/d, c and / being the capacitance and inductance 
per unit length respectively. A typical value is v = 1.2 • 10 8 m/s [40J. Under these conditions 
the Hamiltonian, H = H + Hj, splits into a free part for the qubits and the field 

H = -h{Vt A a z A + Vt A a z A ) + / dkhu k a\a k (2) 

Z J — oo 

and a point-like interaction between them 

Hi = - E djV{xj)o*j (3) 

J=A,B 

Here x j are the fixed positions of the atoms, and dj a x j comes from a dimensional reduction 
of the matter- radiation interaction hamiltonian with two-level atoms and the electromagnetic 
field. We will consider the following initial state 

|^(0)) = M>, (4) 

where only qubit A is excited, while B and the field remain in their ground and vacuum states, 
respectively. We use the formalism of perturbation theory up to the second order and beyond 
Rotating Wave Approximation j8]| and trace over the field degrees of freedom to obtain the 
corresponding two-qubit reduced density matrix p x evaluated at t. In the interaction picture 
with respect to the free Hamiltonian H , the system evolves during a lapse of time t into the 
state 

\m)=T[e- i fi* fExi * )/n ]\e9)®\0), (5) 
T being the time ordering operator. Up to second order in perturbation theory the final state 
can be written as 

m)) = [(1 + A)\eg)+X\ge)}®\0) + 
{U A \gg) + V B \ee))®\l) 

+ (F\eg)+G\ge))®\2) + 0(d 3 ). (6) 

The coefficients for the vacuum, single-photon, and two-photon states, are computed using 
the action (a = A, B) 

St = - l T f V™' (e a \da x a \g a )V(x a ,t')dt' = (7) 
n Jo 

among different photon number states |n),n = 0,1,2..., being \n) {n\ = 
J dk\.... J dk n \ ki...k n ) (k\...k n \ and \k) = a\. |0). Among the various terms present here, 
the only one containing an effective coupling between A and B is 

X = (0\T(S+S A )\0). (8) 
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This includes photon exchange only inside the light cone, vt > r, and vacuum fluctuations 
for all values of t and r, being r = x B — %a the distance between the qubits. The remaining 
terms are 

1 



A =- (0\T(S+S2+SsSi)\0) 



(9) 



;ii^io), 

(2| T(S+S A + S B S+) |0), G = (2| T(StSD |0> . 



U A =(1\SX\0),V B 

fJ- 

2 

Here, A describes intra-qubit radiative corrections, while Ua,Vb,F and G correspond to 
single -photon emission events by one or more qubits. The coefficients in Eq. ([6]) will be 
computed analytically as a function of two dimensionless parameters, £ and K. The first one, 
£ = vt/r, is a dimensionless time variable; the time £ = 1 corresponds to the light-cone, 
which separates two different spacetime regions, before and after photons can be exchanged. 
The second parameter is a dimensionless coupling strength 



K 



Ad 2 N 



(10) 



Note that the qubit-line coupling g = dy/WU/h corresponds to the qubit-cavity coupling 
that appears by taking the same transmission line and cutting it in order to have a length 
L = A (thus creating a resonator). This formulation has the advantage of being valid both for 
inductive and capacitive coupling, the details being hidden in the actual expressions for d and 
N. Tracing over the states of the field, we arrive at the following reduced density matrix 

/ Pu p M \ 

P22 P23 



Px 



(ID 



P23 P33 
V Pu o p 44 / 

representing the two-qubit state in the basis formed by |ee) , \eg) , \ge) , and \gg) . The 
coefficients with the leading order of neglected contributions are 



IVH + O^), P 22 



Pu 

P33 

P14 = U* A V B + 0(d 4 ) 

p 23 = X* + 0(d i ), 
and the state is normalized, c = J2iP 



1 + 2Re(A) + 0(d 4 ) 
X\ 2 + \G\ 2 + 0(d 6 ), P 44 = \u\\ + 0(d 4 ) 



(015+5+10) + 0(d*) 



(12) 



3. Dynamics of Correlations 

In this section we report our results regarding the dynamics of correlations between the two 
(artificial) atoms. We investigate the time evolution of the square root of geometric quantum 
discord yD [|22l . of the entanglement as measured by the negativity iV [27 J, and of the 
maximum connected correlation function C [25] in the state px- We have chosen these three 
specific types of correlations (whose definitions and properties are reported below) for four 
reasons: 
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(i) to best relate our results with the ones reported in [|8l by using a different measure of 
entanglement; 

(ii) to have a more complete description of the time evolution of general quantum 
correlations; 

(iii) to compare quantum correlations with correlations having also a classical nature; 

(iv) to issue a comprehensive comparative analysis among compatible measures of different 
types of correlations. 

The meaning of the latter point can be clarified in connection with our choice of using \J~D. 
One might in fact question why we are comparing different powers of geometric discord and 
entanglement. The reason for this is that we want to be consistent in the order of expansion 
of the perturbative analysis we have performed. A good test to check whether this is true is 
provided by a hierarchy-type relationship that exists between the three chosen quantities for 
arbitrary states p of two qubits, namely 

C{p) > VD(p) > N(p) . (13) 

The rightmost inequality in ( [T3] ) was proven analytically in lfT6ll . while the leftmost one has 
been verified numerically in ||28ll . Notice that for pure two-qubit states both inequalities are 



saturated and ( 13 ) becomes a chain of equalities. 



In our analysis, we have found no violation of the hierarchy (13) for any range of 
the relevant physical parameters characterizing the states px, which serves as a validating 
indication that our results are consistent up to the second order. As a general remark, we 
can conclusively state that all correlations whose time evolution we have looked at, start 
increasing before the time at which the two atoms become causally connected. However, 
the rate of increase changes a lot from one type of correlation to another. 

Let us now introduce the measures of correlations of interest and then discuss the main 
aspects of their dynamical behavior. To this end, the Bloch representation of generic two-qubit 
states p will be useful lfi~9l . Namely, 
1 3 

p = - Rij°i ® °j (14) 

4 i,j=0 



1 

4 



Ii ® h + X] x i a i ® h + Vjh ®<?j + UjCTi <g> Oj 

i=l j=l i,j=l 

where Rij = Tr[p(cr; <g> <x,)], er = I, (i = 1, 2, 3) are the Pauli operators; x = {xi} and 
V = {Hi} represent the three-dimensional Bloch column vectors associated to the qubits A 
and B, respectively; and are the elements of the 3x3 correlation matrix T. 

3.1. Geometric discord 

The geometric measure of quantum discord was first introduced in ll22l and further 
investigations for two-qubit systems have been reported in [|29l [T6l [301 . Given a general 
bipartite 6a® ds quantum state p, the (normalized) geometric discord is defined as 

D(p) = -r^—, min \\p - x\\l (15) 
d.A — 1 xeo 



Dynamics of atom-atom correlations in the Fermi problem 



6 



where x is a so-called classical-quantum state belonging to the set of zero-discord states 
X — 51iPi\i)(i\ <%> QiB and \\P — Q\\l = Tr(P — Q) 2 is the squared Hilbert-Schmidt 
distance between a pair of operators P, Q. We can look at geometric discord as the minimum 
disturbance that would be induced in the system after a projective measurement on one 
of the two parties (say A in the above definition) [|29ll . It is important to remark that its 
value is dependent on the choice of the party to be measured [31 J. Although in principle 
this expression can be very complicated to evaluate explicitly as it involves a minimization 
problem over the set of zero-discord states, an analytical formula exists for the general two- 
qubit case |22l|29l[30]|. In terms of the Bloch picture [Eq. ([T4])], one has 

D{p) = 2Tx[S] - 2A max (S) , 

where A max stands for 'maximum eigenvalue', and the matrix S is defined as S = \{xx T + 
TT T ). The (square root of) geometric discord of px is then 



'D(px) = ^[Re(U* A V B )} 2 + \X\ 2 . (16) 

The above formula is correct up to the second order and it has an immediate physical 
interpretation. The two terms in Eq.([T6|) come indeed from first and second order 
contributions to the time evolution of the state. In particular the X term accounts for 
photon-exchange between the two atoms and carries all the information available about causal 
propagation and atom state dressing. Interestingly, the processes which contribute to non-zero 
quantum discord are and 1 -photon processes and even though the (square root of) geometric 
quantum discord has a continuous evolution starting from t = 0, it is well sensitive to light- 
cone crossing, showing a peak at t = r/c. 



3.2. Negativity 

Negativity ll27l is a well-known and easily computable measure of entanglement for bipartite 
systems which is based on the positivity of the partial transposition (PPT) criterion ll35l . Given 
a general d® d quantum bipartite state p, the (normalized) negativity is defined as 

N(p) = j^-\\p Ta -Iab\\i (17) 

where the T x refers to the partial transposition operation with respect to the x party (x = 
A, B), Iab is the identity operator in the composed Hilbert space "Ha ®%b and ||M||i = 
Tr|M| = J2i \ m i\ is the trace norm for a matrix M with eigenvalues {wij}. As in the case 
of geometric discord we can easily compute the negativity of px up to the second order in 
time-dependent perturbation theory and we obtain the following expression 

N(px) = max jo, ^I^aI 2 - \V B \ 2 ) 2 + W " Wa\ 2 ~ |^| 2 } • (18) 

The three physical processes contributing to entanglement are exactly the same as for 
geometric discord. However, in this case there is a time-dependent condition for entanglement 
to start increasing. Indeed it is easy to check that as long as the following condition is fulfilled 

IXI 2 
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entanglement will be zero. Intuitively, this means that for entanglement to be non-zero, 
second-order processes must dominate over first-order ones. It is worth noticing that such 
a kind of constraint is absent in the case of the (square root of) geometric discord, which 
amounts simply to the sum of two positive and continuous functions. 

3.3. Maximum connected correlation function 

Geometric discord and entanglement are quantities which are strictly connected to the 
quantum character of a system and, indeed, they miss a classical analogue. In this respecr, they 
are key quantities when it comes to understanding the interplay between the foundations of 
quantum mechanics and micro-causality, one of the postulates of relativity theory. However, 
one might also be interested in correlations arising from observable quantities such as, for 
instance, angular momenta and in general spin-like operators. We thus study the atom-atom 
(classical) connected correlation functions Il26ll . to reveal the highest level of sensitivity at 
which the Bloch vectors of the two atoms perceive each other outside and inside the light 
cone. Given a bipartite state p of a pair of two-level quantum systems, we define the maximum 
connected correlation function C(p) as follows 

C(p) = max {((a ■ n) A ® (a ■ h') B ) p - ((a ■ h) A ) p ({& ■ h') B ) p } , (20) 

n,n 

where a is the three-component Pauli-operator vector and (a ■ n) is the projection of such a 
spin vector along the direction pointed by n. For generic two-qubit states p decomposed in 



Bloch form as in Eq. ( |T4| ), the maximum in Eq. ( p0[ ) can be computed in closed form and 
reads [HHHl 

where W = T — xy 1 '. We have computed C{p) for the state p x and obtained an exact 
expression up to the second order, 

C(p x ) = max {\U A \ 2 + \V B \ 2 + 2Re(A), 2(\X\ + \L\)} (21) 

where L = U A V B . Once again, in this case only and 1 -photon processes contribute to the 
above correlation function. Moreover, with C(px) being dependent on X, it shows sensitivity 
to the light-cone crossing. 

3.4. Results and Discussion 

In the following we analyze the time evolution of the above correlations and compare their 
behavior qualitatively and quantitatively. We remark that all of the above quantities depend 
on the three terms U A , V B and X, while the maximum connected correlation function displays 



an A-dependence as well. Fig 1(a) shows the behavior of the square root of geometric discord 
^Jd(px) as a function of the re-scaled time £ = rt/v for different choices of the atom- field 
coupling constant, spanning from a weak to a strong coupling regime, and for a fixed distance 
r = vn / AVL between the qubits. 

The first feature we notice, which is perhaps the most interesting one, is the relatively 
slow but continuous increase that the (square root of) geometric discord shows prior to the 
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Figure 1. (Color online) Time-evolution of (a) the square root of geometric quantum discord 
J D(px), (b) the negativity N(px) and (c) the maximum connected correlation function 
C(px), for r = vir/Afl and three different choices of the coupling strength: Z = 50 (dotted 
green), Z = 200 (dashed red), Z = 400 (continuous blue). 



light-cone crossing. A similar behavior had been previously found by one of the present 
authors in (8l when studying the time evolution of entanglement measured by the concurrence 
O in the same model. However, unlike concurrence which started increasing within a very 
short time interval just before £ = 1, the (square root of) geometric discord reaches finite 
values for a much longer time interval inside the space-like region Js = {0<£<1}. If 
we recall the interpretation of geometric discord mentioned above, we may argue that a one- 
party measurement performed at any time inside Js will always induce a disturbance on the 
composite system. 

Secondly, we observe a peak at £ = 1 which is independent of the interaction regime. 
The height of such a peak, and more generally the global magnitude of the (square root of) 
geometric discord, increases as we increase the coupling. These latter features are easily 
understood by looking at Eq.( 16). As we said above the (square root of) geometric discord 
is the sum of a first order term, which does not carry any causality-related information, and a 
second order term, which instead does carry that kind of information. Hence, the stronger the 
interaction is, the bigger this second order term becomes. 

In Fig 1(b) we show the time evolution of the negativity for the same three choices of 
coupling strengths and the same atom-atom separation. In this case we find essentially the 
same behavior as reported in [8] for concurrence. By comparing Fig 1(a) and 1(b), we may 
conclude that quantum discord is more sensitive to vacuum fluctuations, which are responsible 
for creating correlations between the two atoms outside the light cone. This behavior is well 
understood again when one looks at Eq.(16). The proportionality to X, which is a second- 
order 0-photon term, incorporates exactly this kind of trait. 

Fig 1(c) shows the time evolution of the maximum connected correlation function C(px) 
for the same choice of parameters as in the previous plots. In this case we find something very 
interesting and not at all easily predictable. Indeed, like geometric discord, the maximum 
connected correlation function starts increasing significantly inside Js and it shows a peak at 
£ = 1. The maximum connected correlation function C(p) is not a priori a fully quantum 
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Figure 2. (Color online) Comparative plot displaying the maximum connected correlation 
function C (topmost surface, blue online), the square root \JD of the geometric discord 
(middle surface, red online), and the negativity N (bottommost surface, green online), 
calculated for the state px as functions of the dimensionless time £ and of the coupling strength 

Z, for r = ott/40. 



quantity, and it determines how, on average, the Bloch vectors of the two atoms influence 
each other. The present results seem to suggest that this might be the key quantity to measure 
when it comes to an experimental detection of the dynamics of correlations, provided that a 
simultaneous set of optimal measurements on the two qubits can be efficiently performed in 
the laboratory frame. 

It is worth noticing here that the optimal "measurement directions" h, n' are completely 
different in the space-like region J s and on the light cone. Indeed, for £ < 1, the correlation 
between the two Bloch vectors is best highlighted by measuring the effective spin projections 
in the equatorial x — y plane. On the light cone, on the other hand, the best choice is to 
measure o z for both qubits. This appears to be related to the fact that no excitation can reach 
the atom B before £ = 1 and that, as demonstrated in [8], it is only after this time that the 
excited state population of atom B starts depending on the presence of atom A. The vacuum 
fluctuations, thus, are able to correlate essentially transversal observables for £ < 1, while for 
a longitudinal (z — z) correlation, one has to wait the arrival of the light signal. Finally, in 
Fig. |2]reports a visual comparison of all the three indicators of correlations considered in the 
present analysis, as functions of Z and £. As anticipated, there is no violation of the general 
hierarchy ( [T3| ), thus confirming that the perturbative analysis we have performed is consistent 
up to the present expansion order. 
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4. Non-locality 

One may wonder what the above results mean in term of non-locality. In this respect, a 
key quantity to investigate non-local effects in the dynamics of two two-level systems is 
the Bell parameter 11361 . resulting from well-known inequalities that local classical hidden 
variable theories cannot violate. To make a long story short, one identifies a set of joint 
measurements to be performed on the composite system. Then, based on the outcomes of 
such measurements, it is possible to define a statistical parameter B for which a classical 
threshold value Be exists. Whenever the inequality 

B<B C 

is violated, the state of the system under scrutiny is not reproducible by means of local 
classical hidden variable theories. Unfortunately, the other way around is not true: some 
mixed entangled states exist that do not violate any Bell inequality [37J. In the case of two 
two-level system, it is well known that Be = 2, whereas the maximum violation allowed 
by quantum mechanics is given by the so-called Tsirelson bound, £> max = 2\/2, which is 
saturated by maximally entangled states. Hence, if for a given bipartite quantum state p we 
find 2 < B(p) < G\/g, such a state is not classically reproducible. We have considered two 
different Bell parameters in the present analysis: the conventional CHSH one [38] and its 
optimized version for "X"-shaped states such as p x > given in [|39l . The former, for the states 



in Eq. (11), reads as follows (up to the second perturbative order) 

Bchsh(px) = -\/2(pn + P44 - P22 - P33 + 2Rep 2 3 + 2Repi 4 ) , (22) 
whereas the latter is 11391 

Bopt(px) = 1\fui + max[w 2 ,u 3 ] , (23) 

where 

ui = 4(|pi 4 | + IP23I) 2 «3 = 4(|pi 4 | - IP23I) 2 

U 2 = (Pll + P44 - P22 - P33) 2 • 

The above quantities correspond to two different choices of the Bell parameters, that is, 
different choices of the angles along which we project the effective spin operators in a joint- 
measurement experiment. The main difference between them is that Bopt is optimal in the 
sense that it maximizes the violation of the related Bell inequality, whenever such a violation 



is present. We report in Fig 3(a) the time evolution of the the Bell parameter Bqhsh- It 
is clear from these plots that, in order to observe a violation of the Bell inequality, a very 
strong coupling is required (Z ps 1000). However, such a violation is witnessed only in the 



surroundings of the light cone crossing £ = 1. Fig 3(b) shows, instead, the time evolution 
of Bopt- Some qualitative and quantitative differences between these quantities, especially 
for £ > 1 and very strong coupling, are present. First of all we notice that, as expected, the 
optimized Bell parameter Bopt is greater than Bqhsh for all couplings and times. Secondly, 
the latter is clearly more sensitive to the strong-coupling regime. This result makes perfectly 
sense if we look at the entanglement dynamics as a function of the coupling constant: the 
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Figure 3. (Color online) (a) Bell-CHSH parameter Bchsh and (b) its optimized version 
Bqpj, plotted for Z = 100 (dotted green), Z = 400 (dashed red) and Z = 800 (continuous 
blue), as a function of £. The straight line at B = 2 gives the limit for a local realistic 
description. As one can see, in both cases, a very strong coupling is needed in order to 
have a clear violation of the inequality, occurring, however, only very close to the light cone 
crossing point, (c) Comparison between Bchsh (red dashed) and Bqp^ (blue dotted) for 
strong coupling Z = 1000. In all these plots r = vtt/AQ. 



bigger Z is, the more entanglement is present in the system, pushing up the Bell parameters' 
value. These last two features are best enlightened in the third plot |3(c)| where the dynamics 
of the two Bell parameters we considered, is compared in the case of very strong coupling 
Z = 1000. 



5. Experimental implementation 

Now we will focus on the particular experimental set-up that we propose to test the results 
above. Clearly, we need the ability of fast on-off switching of the qubit-field interaction. 
We should be able to prepare the initial state Eq. ([4]) and let the interaction be active only 
during a finite time. We can conceive several circuit QED schemes in order to achieve these 
goals. As commented before our general formalism can accommodate several architectures; in 
particular, it is valid both for inductive and capacitive couplings. We can choose for instance 
a pair of three-junction flux qubits galvanically coupled to the center conductor of an open 
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transmission line. In Ref. [7] it is shown how the desired initial state can be prepared with high 
fidelity by varying an external magnetic flux adiabatically for qubit B and non-adiabatically 
for qubit A. After that, the interaction has to be switched on and kept constant during a given 
time interval. In Ref. IPTlII several modifications of the three-junction scheme are proposed in 
order to achieve couplings tunable in strength up to the ultra-strong regime. In particular, a 
specific set-up featuring an intermediate superconducting loop has been described in details 
in Ref. [|42l . where an effective interaction Hamiltonian between a qubit and a transmission 
line has been derived, that reads 

Hj(xcos(f)a x V(x), (24) 

/ being related with an external magnetic flux $ threading the SQUID: / = 27r$/$ (where 
$o = h/2e is the flux unit). By a suitable modulation of the flux $, thus, the interaction can 
be activated and switched off at will. State of the art circuit-QED technology allows variations 
of the magnetic flux at frequencies of about 10 GHz [|40l and larger values are expected for the 
future 11431 . As a result, switching times of the order of 0.1 ns and even shorter can be safely 
considered. Taking the values that we have considered in our plots, and typical circuit-QED 
parameters, such as v = 1.2 ■ 10 8 m/s and = 10 9 m/s, the point £ = 1 is equivalent to an 
interaction time t ~ 1 ns. Thus, the region around £ = 1 is well within experimental reach. In 
particular the strongest value that we considered for the adimensional coupling Z is equivalent 
to g/Cl ~ 0.3, which is quite similar to the ones in cutting-edge experiments investigating the 
ultra-strong coupling regime [|44ll451 . 

Once the interaction is switched off, quantum state tomography [46J may be performed 
in order to quantify the degree of correlations, using for instance the magnitudes that we have 
considered in this work. The dynamics is effectively frozen and the system remains in the 
state px{t), so measurements can take as much time as required. In particular, it would be 
interesting to run the experiment for different interaction times inside and outside the light 
cone, in order to test both the peak at £ = 1 and the correlations for £ < 1. 

6. Concluding Remarks 

To summarize, using second order perturbation theory we have discussed the dynamics of 
quantum correlations in the Fermi problem, which can be experimentally tested in a one- 
dimensional set-up involving two artificial atoms coupled to the electromagnetic field of an 
open ended transmission line. We have compared the time behavior of the entanglement, as 
measured by the negativity, and of more general quantum correlations such as the (square 
root of) geometric discord, and the maximum connected correlation function. All of these 
correlations display a peak at the light cone crossing point, £ = 1, which corresponds 
to the time at which a signal from atom A arrives at B. The geometric discord and 
the connected correlation, however, have a substantially non-zero value also in the space- 
like region. This is due to the fact that electromagnetic vacuum fluctuations can induce 
(transverse) correlations that are signalled by these functions. As the light-cone is crossed, 
these correlations change their character and become longitudinal as a result of the fact that 



Dynamics of atom-atom correlations in the Fermi problem 



13 



the excited state population of the second atom starts to depend on the presence of the first one. 
We have briefly investigated non-local effects in this model as encoded in possible violations 
of Bell inequalities. We have found that a violation can occur in the neighborhood of the light- 
cone crossing and only for strong couplings between the atoms and the propagating field. We 
believe that both geometric discord and maximum connected correlations, which are found 
to be sensitive to causal propagation, are suitable candidates for understanding and testing 
experimentally the role of micro-causality in the dynamics of quantum correlations. 
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